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Abstract 

The Goldstein-Ruban theory has been extended 
within the quasi-parallel framework of ZavoEskii et al. 
to study the acoustic receptivity of compressible bound- 
ary layers. We consider the receptivity produced in a 
region of localized, small-amplitude variation in the sur- 
face temperature and compare it with the receptivity that 
is induced through a similar mechanism by a variation 
in the suction velocity at the surface. The orientation 
of the acoustic wave can have a significant impact on 
the receptivity process, with the maximum receptivity 
at a given sound-pressure level produced by upstream- 
oriented acoustic waves. At sufficiently low Mach num- 
bers, the variation of receptivity with the acoustic-wave 
orientation can be predicted analytically and is the same 
for both surface suction and surface heating. However, as 
a result of the acoustic refraction across the mean bound- 
ary layer, the above dependence can become complex and 
dependent upon the type of surface nonuniformity. The 
results also suggest that the receptivity caused by tem- 
perature nonuniformities may turn out to be more signif- 
icant than that produced by the mean-flow perturbations 
associated with strip suction. 

1. Introduction 

The central issue in regard to the generation of insta- 
bility waves in low-speed boundary layers was the iden- 
tification of the wavelength-reduction mechanisms that 
enable the long-wavelength free-stream disturbances to 
excite short-wavelength instability waves (Reshotko 1976, 
Goldstein 1983). The pioneering work of Goldstein (1983, 
1985), as well as that of ZavoEskii et al. (1983), Ruban 
(1985), and Goldstein et al. (1987), showed that the 
key to this wavelength-reduction process is the occur- 
rence of sufficiently shorter length scales in the mean flow, 
which can tune the free-stream disturbance to the length 


scale of the boundary-layer eigensolutions. These short- 
scale variations may occur either naturally (because of 
the rapid development of a boundary-layer flow in the 
vicinity of the leading edge (Goldstein 1983)), or may 
be induced externally (by a short-scale variation in the 
surface boundary conditions). Examples in the second 
category include irregularities in the shape of the airfoil 
surface (ZavoEskii et al. 1983, Goldstein 1985, Ruban 
1985) and nonuniformities in the surface distributions of 
suction velocity and surface admittance (Kerschen and 
Choudhari 1985). The theories of Goldstein, ZavoEskii et 
al . , and Ruban have subsequently been applied to a vari- 
ety of problems, and the reader is referred to the reviews 
by Goldstein and Hultgren (1989) and Kerschen (1989), 
as well a s to the relevant papers in Reda et al. (1991). 
Many significant contributions in the area of boundary- 
layer receptivity have been reported in the Russian litera- 
ture as well, and some of those papers have been reviewed 
by Nishiokaand Morkovin (1986) and Kozlov and Ryzhov 
(1990). 

Much of the above research was related to the acous- 
tic generation of instability waves and focused on the 
incompressible flow regime, in which the wavelength 
of acoustic disturbances is effectively infinite. The 
acoustic wavelength shortens progressively (in compar- 
ison with the convective length scale at the same fre- 
quency) as the Mach number is increased. However, 
it remains asymptotically larger than the wavelength 
of the viscous-inviscid interactive Tollmien-Schlichting 
mode (TS) (Ryzhov and Zhuk 1986, Smith 1989) which 
dominates the primary instability of two-dimensional 
boundary layers at subsonic and moderately supersonic 
speeds. Consequently, the wavelength reduction mecha- 
nisms previously identified by Goldstein (1983, 1985), and 
Ruban (1985) are also directly relevant to the generation 
of the viscous-inviscid interactive modes in compressible 
boundary layers. 



At high supersonic and hypersonic Mach numbers, 
the dominant instabilities of a boundary-layer flow are 
predominantly inviscid, and their phase speeds (or wave- 
lengths) are comparable to those of the acoustic distur- 
bances in the free stream. However, because these invis- 
cid instabilities (i.e., Rayleigh modes) are primarily sub- 
sonic with respect to the local free stream, their phase 
speeds do not match with the acoustic propagation speed 
(except in a limiting sense as described in the following 
paragraph). Therefore, a spatial tuning of the unsteady 
motion (forced by the free-stream disturbances) is still 
necessary for the excitation of the inviscid instabilities. 
This tuning can be provided by the short-scale surface 
nonuniformities in much the same manner as that in the 
incompressible boundary layers. However, the asymp- 
totic scaling involved is completely different for the in- 
viscid modes of supersonic boundary layers. Accordingly, 
the details of the energy transfer process can be signifi- 
cantly different in the two speed regimes, as discussed by 
Choudhari and Streett (1990). This work describes the 
extension of the Goldstein- Ruban theory to high-speed 
boundary layers and presents numerical results for both 
viscous-inviscid interactive modes and predominantly in- 
viscid Rayleigh instabilities. However, because the main 
purpose of that paper was to provide an overview of the 
receptivity of high-speed boundary layers, only a few sets 
of results were presented. 

The limiting case mentioned in the preceding para- 
graph involves a resonance in the leading-edge region of 
a flat-plate boundary layer between the inviscid insta- 
bility modes and the acoustic free-stream disturbances 
that propagate parallel to the plate (Gaponov 1985). The 
strong excitation of instability waves by such “grazing” 
acoustic disturbances was recently explained by Fedorov 
and Khokhlov (1991) with a theory that is the supersonic 
analog of Goldstein’s (1983) seminal work on TS-mode 
generation near the leading edge of a flat-plate airfoil. 
Because of the above-mentioned resonance feature, the 
study of leading-edge receptivity in a highly supersonic 
boundary layer reduces to a study of the sound-generation 
mechanisms that excite the grazing acoustic disturbance 
via a scattering process near the leading edge. (See the 
more recent work of Fedorov and Khokhlov (1993) in this 
context.) 

In addition to the receptivity produced near the lead- 
ing edge and the receptivity caused by scattering near 
surface nonuniformities, other possible mechanisms ex- 


ist through which a TS or a Rayleigh instability can be 
excited in a supersonic boundary layer. These mecha- 
nisms include a resonant forcing through the quadratic 
interaction between a pair of acoustic and convected free- 
stream disturbances, and a sustained generation of three- 
dimensional TS waves through the interaction of slow, 
near-planar acoustic disturbances and long-wavelength 
Gortler vortices. However, these mechanisms are some- 
what specialized, and one may expect the overall receptiv- 
ity to be dominated by the leading-edge mechanism and 
the surface nonuniformities. The relative importance of 
the two dominant routes for receptivity will depend cru- 
cially upon the Mach number of the flow, as well as on 
the leading-edge geometry and the amplitude and spatial 
distribution of the surface nonuniformities. This particu- 
lar issue is beyond the scope of this paper and, therefore, 
will be addressed in a future article. 

In this paper, we will examine in further detail the 
generation of viscous and inviscid instabilities via the in- 
teraction of a time-harmonic, acoustic disturbance with 
short-scale small-amplitude variations in the surface dis- 
tributions of suction velocity and temperature. The ef- 
fectiveness of strip suction and strip heating in stabilizing 
gas boundary layers at subsonic speeds was demonstrated 
theoretically by Reed and Nayfeh (1985) and Masad and 
Nayfeh (1992), respectively. However, the mean-flow dis- 
turbance induced by these devices can also act as a cat- 
alyst in the receptivity process, as described above. The 
receptivity due to both localized and nonlocalized varia- 
tions in the suction-velocity distribution was studied in 
the context of low-speed flows by Kerschen and Choud- 
hari (1985) and Choudhari (1992), respectively; Choud- 
hari and Streett (1990) investigated the localized recep- 
tivity in high-speed flows. The use of surface heating or 
cooling for laminar flow control (LFC) is primarily rele- 
vant to the case of high-speed boundary layers, and the 
associated receptivity mechanism does not appear to have 
been studied thus far. 

As shown by Kerschen and Choudhari (1985), two 
distinct scattering mechanisms are responsible for the 
acoustic receptivity induced by a porous surface that is 
typically used for suction-based LFC. These two mech- 
anisms are (1) scattering by the short-scale mean-flow 
variations produced by a specified surface-suction distri- 
bution and (2) scattering by variations in the acoustic 
admittance of the porous surface. When the amplitudes 
of both variations are sufficiently small, the receptivity 


2 



produced via each mechanism can be examined indepen- 
dently of the other mechanism. In this paper, we focus 
on the receptivity produced in the zero-admittance limit 
(i.e., on the receptivity that occurs through the mean- 
flow variations only). A comparison will be made with 
the analogous receptivity mechanism in the case of heat- 
transfer strips at the surface. Receptivity that occurs via 
surface admittance variations is studied in a separate pa- 
per, wherein we also examine the effect of a nonzero sur- 
face admittance on the stability of a compressible bound- 
ary layer. 

Generally, in this paper we focus on the localized 
receptivity process associated with a single suction, heat- 
ing, or cooling strip. However, receptivity produced in 
noncompact LFC configurations that involve a series of 
such strips will also be considered briefly. To predict the 
outcome of the localized receptivity process, we utilize a 
nonasymptotic extension of the Goldstein- Ruban theory, 
which was summarized by Choudhari and Streett (1990). 
A similar approach was used earlier by Choudhari and 
Streett (1991) to predict the suction-induced receptiv- 
ity in low-speed boundary layers. The assumption of a 
small-amplitude surface disturbance may appear to be 
somewhat restrictive in terms of applying the theory in 
practice. However, this assumption simplifies the analysis 
to a great extent; moreover, it enables fairly general con- 
clusions to be drawn in regard to the receptivity produced 
by different possible combinations of free-stream and sur- 
face disturbances, irrespective of the particular geometry 
involved. (See Goldstein 1985.) A brief description of the 
above nonasymptotic approach applied to compressible 
boundary layers is given in section 2 below. The numer- 
ical results, which are obtained for a flat-plate boundary 
layer at various Mach numbers, are presented in section 
3. Conclusions are given in section 4. 

2. Description of the Analysis Involved 

Consider the uniform two-dimensional steady super- 

sonic flow of a perfect gas past an infinitely thin and 

nominally flat airfoil that has a localized nonuniformity 
in the surface distributions of suction velocity and/or the 
temperature at a distance t from the leading edge. The 
velocity of the oncoming free stream is denoted by , 
and the corresponding Mach number is symbolized by 
M. The pressure, temperature, density, absolute viscos- 
ity, and specific heat at constant pressure at the free- 
stream conditions are denoted by , f.i ^ , and 


<7* , respectively; the ratio of the specific heat at con- 

stant pressure to that at constant volume is denoted by 7. 
We now assume that the local distributions of the surface 
suction velocity and the surface temperature are given by 

V:/U^=(^Fi s \X,Z) , (la) 

and 

TJIZ = T W:0 +c$Ft t \X,Z) , (lb) 

respectively. The amplitude parameters and fl/ are 
assumed to be sufficiently small, so that the mean-flow 
disturbance produced by each type of nonuniformity can 
be treated as a small perturbation to the unperturbed 
mean boundary layer. The quantity T w q in (16) repre- 
sents the local (nondimensional) surface temperature for 
the unperturbed mean flow. The nondimensional func- 
tions Fw\X } Z) and F^\x, Z) denote the normalized 
distributions of suction velocity and surface temperature 
in terms of the local coordinates 

X = (x* - t)/6* and Z = z*/6* (1 e) 

along the streamwise and spanwise directions, respec- 
tively, where <5* is the local displacement thickness as- 
sociated with the unperturbed mean boundary layer. 

The unsteady free-stream disturbance incident upon 
the surface irregularity is assumed to be a time-harmonic 
acoustic wave with the frequency u)* (comparable to the 
range of locally unstable frequencies), the streamwise 
wavenumber a* c , and the spanwise wavenumber /3* c . We 
use the symbol c ac to denote the nondimensional ampli- 
tude of this wave, which is assumed to have been defined 
in terms of the amplitude of the acoustic fluctuation in 
an appropriate flow variable and normalized by a suitable 
quantity associated with the mean flow. For now, we will 
leave open the choice of the flow quantities involved in 
defining ( ac and will simply assume that the latter pa- 
rameter is sufficiently small, so that the unsteady part of 
the motion can always be treated as a small perturbation 
to the local mean flow throughout the region of interest. 

As first shown by Goldstein (1985) and Ruban 
(1985), the problem of studying the interaction of the 
acoustic wave with the steady nonuniformity on the air- 
foil surface is simplified considerably after the the small- 
disturbance approximation (eL 2 * * 5 \^ << 1, <ac << I )is 
invoked and the local nature of the interaction pro- 
cess ( 6 * << £*) is exploited. Specifically, if we use 
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UZoiPloVZo 2 * an ^ TZo to nondimensionalize the local ve- 
locities, pressure, and temperature, respectively, and 
denote the corresponding nondimensional quantities by 
dropping the superscript * on them, then we can expand 
the array of dependent variables Q = ((/, V } W, T, P) for 
each type of surface nonuniformity in terms of the two- 
parameter perturbation series 

Q W(X,Y,Z,t) = Qo (Y) + tV>QHXX t Y t Z) 

+ fa C Qac(V')e*'b*'' v+ ^ z - w '] 

+ % ) f a cQ$„(X,Y,Z)e- iut 

+0(^ )2 ,4). (2a) 

where the superscript j indicates whether the surface 
nonuniformity involves a short-scale variation in the suc- 
tion velocity (j = s ) or in the surface temperature ( j = t). 
The wall-normal coordinate Y has been scaled with re- 
spect to 6*, whereas the frequency u> and time t have 
been normalized by U^/6* and its inverse, respectively. 
Implicit in the form of (2a) is the assumption that the 
unperturbed boundary-layer flow is convectively unsta- 
ble (Huerre and Monkewitz, 1990) at the location of in- 
terest. This assumption is central to most problems of 
receptivity, as discussed by Goldstein (1985). 

The nonzero elements of the zeroth-order term 

Qo(V') = (U 0 (Y), 0, 0, T 0 (Y), ft—) (26) 

denote the quasi-parallel approximations to the corre- 
sponding base-flow quantities at x* = t* . The first-order 
perturbations due to the wall nonuniformity and the free- 
strearn acoustic wave are indicated by the subscripts w 
and ac, respectively. We can easily show that the mean- 
flow disturbance due to the surface nonuniformity 
satisfies the steady version of the parallel-flow disturbance 
equations (Lekoudis et al. 1976), subject to inhomoge- 
neous boundary conditions at the surface that are ap- 
propriate for the type of nonuniformity involved. It is 
necessary to point out that the surface-temperature vari- 
ation induces a change in the local mean flow through 
changes in both the density and the viscosity of the flow. 
In this paper, we will not make any attempt to assess 
the separate contributions to the mean-flow perturbation 
through the above two channels; however, we expect that 
for M = 0(1) the dominant contribution will be provided 
by the changes in the density of the flow. The signature 


Q ac of the acoustic disturbance within the mean bound- 
ary layer satisfies the unsteady parallel-flow disturbance 
equations, subject to outer boundary conditions appro- 
priate for the specified form of the incoming disturbance 
(Mack 1984). 

Neither of the two first-order perturbations simul- 
taneously possesses the combination of length and time 
scales relevant to the local instability motion. However, 
this required combination is produced by the scattering 
of the acoustic-signature field by the mean-flow distur- 
bance. The leading-order solution for this scattered field 
is given by the 0(e«?^ac) term in (2a) (i.e., the term 
that contains Qw,ac)j which represents the bilinear in- 
teraction between the first-order perturbations and 
Q ac . After a Fourier transform is taken along the stream- 
wise (X — ► a) and spanwise ( Z — ► ft) directions, 
satisfies an inhomogeneous form of the linear stability 
equations 

L NS (a,dY,(i,uj,R 6 .)QW ac = Njvs ;a (Q«>,Q a(! ), (3) 

in which the overbar denotes the Fourier transformed 
quantities. 

The operator L on the left-hand side of (3) de- 
notes the matrix form of the linearized continuity and 
Navier-Stokes equations in the transform space (cv, ft), 
and ~ is the Reynolds number based 

on the local displacement thickness of the unperturbed 
boundary layer. The expression for Lyvs can be found in 
different equivalent forms in the literature on the quasi- 
parallel stability of compressible shear flows. The forcing 
vector N N5i 2 on the right-hand side represents the bilin- 
ear interaction (indicated by the subscript 2) between its 
vector arguments and Q ac . It arises from the con- 
vective and viscous nonlinearities in the exact governing 
equations. The expressions for the elements of Nyvs ; 2 are 
provided in the appendix. Even though we are interested 
only in the receptivity of a two-dimensional base flow in 
this paper, the expressions given in the appendix are gen- 
eral and may be applied to arbitrary three-dimensional 
base flows and other types of free-stream disturbances. 
In fact, with a suitable modification of the terms that in- 
volve w, the same expressions may be utilized in a study 
of resonant triad interactions among disturbaces in com- 
pressible shear flows. As pointed out by Choudhari and 
Streett (1990), similar interactions that involve an ap- 
propriate combination of an acoustic and a converted 
free-stream disturbance could also lead to receptivity in 
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high-speed boundary layers. However, the receptivity via 
this latter mechanism is expected to be relatively weak 
in view of the quadratic dependence of receptivity on the 
amplitude level of the free-stream disturbances, which is 
typically very small. 

Depending on the geometry F^\X,Z) of a three- 
dimensional surface nonuniformity, the instability motion 
excited in its vicinity can have a considerably complex 
spatial structure. However, because the interaction that 
leads to this receptivity process is bilinear, the instabil- 
ity waves that are produced via each spanwise Fourier 
component of F^\X,Z) can be analyzed independently. 
Therefore, we will confine our attention for the remain- 
der of this work to those surface disturbances that are 
periodic in the Z direction with a (single) wave number 
ft = ft w . That is, we will focus on cases where 

F {i \X,Z) = T {i \X) e'P™ 2 . (4a) 

The numerical results to be presented in section 3 cor- 
respond exclusively to two-dimensional surface distur- 
bances (i.e., to the special case ft w = 0). In view of (4a), 
the unsteady scattered field Q &?ac has the form 

Q$ac(X> y > Z) = Y ) (46) 

The part of Q^'' ae that corresponds to the instability 
mode of interest can be found as the residue contribu- 
tion to the inverse Fourier integral 

1 f°° 

QW*e(X, Y )=-*r Y)e' aZ da (5a) 

from an appropriate singularity of the transform solution 
Y) in the streamwise wave-number (a) space 
(Goldstein 1985, Ruban 1985). In general, this singularity 
assumes the form of a first-order pole, and its location is 
given by the spatial dispersion relation 

ot = (*ins(fiins,W, Rf) {Pins = Pac + Pw) ( 56 ) 

for the instability mode under consideration. As a re- 
sult of isolating the instability-wave part of the unsteady 
scattered field, one is able to relate the amplitude of the 
generated wave to the amplitude of the incident acoustic 
wave. Thus, the dimensional perturbation in the stream- 
wise component of the unsteady mass flux that is associ- 
ated with the generated instability wave can be expressed 
in the form (Goldstein 1985, Choudhari and Street! 1990) 


(pUY in ,(X,Y,ZJ) = 

(pVYac.i C pU,pU E po(Y) , (6a) 

where pU* ci is the streamwise component of the mass- 
flux perturbation that is produced by the incoming acous- 
tic disturbance at the airfoil location (but in the absence 
of the airfoil). The quantity E p u(Y ; u>, P w +p aC i Rf) rep- 
resents the local eigenfunction for the pU perturbation, 
which is assumed to have been normalized so that it has a 
maximum magnitude of unity across the boundary layer. 
The “coupling coefficient” function (Tam 1981, Goldstein 
1985) pU {0 w , a ac , /3 ac ,u , R s -) then represents the ra- 
tio of the amplitude of the generated instability wave at 
X = 0 to the amplitude of the incident acoustic wave; 
both amplitudes are measured in terms of the respective 
fluctuations in the streamwise mass flux pU , as indicated 
by the double subscript pU ; pU . 

When the geometry of the surface disturbance is 
more general than the harmonic form in (4a), the to- 
tal unsteady fluctuation associated with the locally ex- 
cited instability motion can be evaluated by integrating 
over the contributions from the different spanwise Fourier 
modes of the surface nonuniformity. (See Choudhari and 
Kerschen (1990) and Choudhari (1993).) Thus, if the 
nonuniformity has a periodic structure in the spanwise 
direction (e.g., that corresponds to a periodic array of 
suction holes), then the above sum will assume the form 
of a discrete Fourier series in the Z coordinate. On the 
other hand, if the nonuniformity is localized in both the 
streamwise and spanwise directions, then this sum will 
become an inverse Fourier integral in the fti ns plane. 

Because of the bilinear nature of the acoustic scat- 
tering process, the coupling coefficient function can be 
written in a form that is similar to that in low-speed 
boundary layers, namely, 

^pU\pV ~ e tu^ F^\ q w) ; p u ( w » ,°L aC) p ac ), (66) 

where a w = a ins - a ac and the cases j = s and j - t 
correspond to the receptivity induced by variations in 
the wall suction velocity and the wall temperature, re- 
spectively. The factor measures the amplitude 

of the Fourier harmonic (ct w ,ft w ) of the surface suction 
or the surface temperature distribution. This particu- 
lar Fourier component of the mean-flow disturbance is 
responsible for “tuning” the acoustic phase to that of 
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the generated instability mode. The efficiency function 

p v ls independent of the local geometry and is com- 
pletely determined by the type of surface nonuniforrnity 
and other flow parameters such as the frequency u and 
wavenumbers c* ac and j3 ac of the incoming acoustic wave, 
and the local Reynolds number R <$• . This type of decou- 
pling between the geometry and efficiency factors, which 
was first pointed out by Goldstein (1985), is common to 
all receptivity mechanisms that involve a scattering of 
small-amplitude free-stream disturbances by weak surface 
inhomogeneities. As the numerical results in the follow- 
ing section demonstrate, the above decoupling enables 
(1) the study of the intrinsic features of each receptivity 
mechanism independent of the geometry involved and (2) 
the direct comparison of the receptivity caused by differ- 
ent types of surface nonuniformities, with the assumption 
that the spatial distributions for all nonuniformities are 
similar to each other (but otherwise arbitrary). 

In (6), we used the perturbation in the streamwise 
mass flux to measure the amplitudes of both the input 
(i.e., acoustic) and the output (i.e., instability) waves. As 
noted by Choudhari and Streett (1991), the knowledge of 

P u [S adequate to determine the outcome of the lo- 
calized receptivity process in any given case. However, 
the trends in the coupling-coefficient magnitude with re- 
spect to variation in the flow parameters can be substan- 
tially different if either or both of the wave amplitudes are 
measured in terms of a different physical quantity, such 
as the associated pressure fluctuation at the surface. The 
selection of the most appropriate physical quantity for 
measuring the amplitudes of the reference and the out- 
put waves is relatively less difficult in the case of incom- 
pressible boundary layers, particularly when these flows 
are two dimensional or axisymmetric. This is because the 
dynamics of the transition process in the above cases is 
dominated by velocity perturbations that are parallel to 
the mean flow. 

The streamwise velocity perturbation is also easy to 
measure in a wind-tunnel experiment by using a hot-wire 
anemometer. Because the hot-wire response is actually a 
combined measure of the perturbations in both the veloc- 
ity and the density at the measurement location, a natu- 
ral extension of the above definition of wave amplitudes 
to compressible boundary layers is to specify the wave 
amplitudes in terms of the respective fluctuations in the 
streamwise mass flux. This choice was made by El-Hady 
(1991), who examined the effect of flow nonparallelism on 


the growth of an instability wave. In studies of receptivity 
in compressible flows, the same choice of wave amplitudes 
is deemed convenient under many circumstances. How- 
ever, situations do occur in which the pressure fluctuation 
produced by the incident acoustic wave is more appropri- 
ately used as the measure of the input-wave amplitude. 
This situation occurs, for instance, when one is study- 
ing the variation of the coupling coefficient with respect 
to the orientation of the incident acoustic wave because 
the pU perturbation associated with an incident wave of 
given sound pressure level (SPL) can vary dramatically 
as the angle of incidence is varied across its permitted 
range. 

3. Numerical Results 

For a parametric study, we choose to examine the 
receptivity of the boundary-layer flow past a semiinfi- 
nite, adiabatic flat plate. The stagnation temperature 
of the oncoming flow is held fixed at 311 K as the Mach 
number is varied. To start, we will consider the gen- 
eration of a two-dimensional instability wave (by two- 
dimensional acoustic disturbances) in a subsonic bound- 
ary layer. Specifically, in section 3.1, we examine the in- 
fluence of the free-stream Mach number on the coupling 
coefficients for receptivity that is produced through vari- 
ations in both surface suction and surface temperature. 
Of particular interest is the comparison of the receptiv- 
ity produced by the two types of LFG devices for a given 
extent of stabilizing influence on the incoming boundary 
layer. The influence of acoustic-wave orientation on the 
coupling coefficients is considered in section 3.2. Therein, 
we illustrate how the acoustic-signature field and the cou- 
pling coefficient in each case can at times be dramatically 
influenced by the refraction of the acoustic motion across 
the mean boundary layer. The cumulative effect of a 
series of suction or heating or cooling strips on the recep- 
tivity process is studied in section 3.3. In subsonic flows 
at sufficiently large Mach numbers, the spanwise length 
scale of the acoustic waves can become sufficiently short, 
whereby a three-dimensional instability wave can be ex- 
cited even when the surface nonuniformity (i.e., suction 
or heating or cooling strip) is purely two dimensional. 
The excitation of three-dimensional instability waves by 
three-dimensional acoustic waves is examined in section 
3.4. Finally, in section 3.5, the receptivity of a supersonic 
boundary layer is considered. 
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3.1 Influence of Mach Number on the Cou- 
pling Coefficient 

In the first set of calculations, we consider the gen- 
eration of two-dimensional instabilities near a suction (or 
heat-transfer) strip that corresponds to a uniform distri- 
bution of suction velocity (or temperature perturbation) 
across its width. The location of the strip is assumed to 
be at R = R l f l 2 “ 955 at all Mach numbers. This choice 
results in a slightly different value of R&* at each Mach 
number. However, the use of Rf (i.e., the Reynolds num- 
ber based on the distance from the leading edge) is more 
convenient for a flat-plate geometry. Moreover, because 
the range of Mach numbers considered here is limited, the 
accompanying variation in Rs* is not large. 

We assume that the orientation of the incident acous- 
tic wave is held fixed at 0 ac = 45° as the Mach number 
is varied. Here, 0 ac denotes the (clockwise) orientation 
of the wave-number vector of the incident acoustic wave 
with respect to the mean-flow direction. Thus, the range 
0 ac £ [0,tt/ 2) corresponds to incident waves that travel 
downstream, whereas 0 ac £ (tt/ 2, 7 r] denotes the range of 
upstream-traveling incident waves. The polar angle of in- 
cidence 0 v ac with respect to the surface normal is related 
to 0 ac via 6? c = |7t/2 - 0 ac |. Note that in a wind-tunnel 
experiment it is perhaps easier to produce an acoustic 
wave that propagates downstream (as well as upstream 
in a subsonic flow) in a direction parallel to the plate. 
However, the choice 6 ac = 45° appears to be more gen- 
eral than 9 ac — 0° for assessing the values of the coupling 
coefficient. Subsequent calculations in which we examine 
the influence of acoustic-wave orientation on receptivity 
will demonstrate that the value of the efficiency function 
for an incident wave that propagates downstream (and 
upstream if the Mach number is small) can be easily re- 
lated to the efficiency function for an incident wave that 
propagates along 9 ac = 45°. 

The geometry factor for the suction or heat-transfer 
strip is given by 

pi) = 2 1/2 sinK,rf/2) (7a) 

ir ] / 2 ct w 

where d is the (nondimensional) width of the strip. We 
have assumed that the amplitude of T^\X) is equal to 
unity for both j = s and j = t. We will focus our atten- 
tion on the excitation of the instability mode that corre- 
sponds to the lower branch (to be indicated by subscript 
lb in this paper) of the neutral stability curve at R = 955. 


In that case, (7a) shows that the geometry factor has its 
maximum magnitude when the width d of the suction or 
heating strip is given by 

7 f . 

d — (^u; ; lb — ^ac; lb ) * 0 

«u/; ib 

To obtain the upper limit on the receptivity produced by 
a single strip, assume that the strip width is equal to the 
value given by (76). Because the wave number «i„ 5; ib 
of the neutral instability wave at R = 955, and the wave 
number a ac> ib of the incident acoustic wave vary with M, 
the total mass or heat flux across the width of the strip 
will also be a function of M with the above choice of strip 
width. This total-flux parameter is expected to determine 
the overall stabilizing influence of a suction (or heating) 
strip upon the boundary layer. As a result, the total 
“gain” from the suction or heating strip will vary to some 
extent with the Mach number. Nevertheless, the overall 
variation in a w -\b is relatively small at the assumed value 
of acoustic orientation («u, ; ib ~ 0.13 at = 0, wffiereas 
cx w - ib « 0.09 at M - 0.99). Therefore, the above choice of 
strip width appears to be reasonable from the standpoint 
of assessing the effect of free-stream Mach number on the 
coupling coefficient C^ s J. pU that is likely to be measured 
in an experiment. 

In figures 1(a) and 1(b), we have plotted the magni- 

C<0 

tudes of the normalized coupling coefficients and 

C (0 

pU ( i( LI in the above case as functions of M . For both 
types of surface disturbances, the coupling-coefficient 
magnitude varies linearly with M as M — ► 0. In other 
words, the compressibility correction to the coupling coef- 
ficient is O(M) in magnitude, which is significantly larger 
than the 0(M 2 ) variation in the same limit in both the 
local stability properties (Ryzhov and Zhuk 1986, Smith 
1989) and the characteristics of the mean-flow pertur- 
bation (Stewartson 1974). The O(M) variation in the 
receptivity properties is mainly caused by a comparable 
variation in the acoustic-signature field as M — ► 0 and 
can be explained in the following manner. 

Goldstein (1985) and Ruban (1985) showed that the 
receptivity near the lower branch location is mainly con- 
centrated in a thin subregion of the mean boundary layer 
that occupies the region close to the surface. Because the 
acoustic wavelength at the frequency and in the range 
of Mach numbers under consideration is asymptotically 
larger than the thickness of the mean boundary layer, the 
acoustic pressure fluctuation varies very little across the 
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boundary layer. Moreover, the zero-normal-flow condi- 
tion at the surface causes the reflected pressure wave to 
be in phase, locally, with the incident wave. As a result, 
the pressure fluctuation near the surface is twice as large 
as the pressure perturbation associated with the incident 
wave alone. The pressure perturbation drives a Stokes 
shear wave in the viscous wall layer, which is scattered 
by the local mean-flow perturbation to produce the in- 
stability wave. The profile of the Stokes wave is indepen- 
dent of the flow Mach number, the acoustic orientation 
9 aC) and the mean boundary-layer properties. As a conse- 
quence, the coupling coefficient is simply proportional to 
the amplitude U s t of the Stokes wave, which differs from 
the streamwise velocity perturbation associated with the 
incident acoustic wave by a factor given by 


U»t 

U a c : t 


2 

1 + M cos 9 a 


(^ac; * — ISA COS 0 a c P a c ; I ) * 


( 8 ) 

The factor 2 in the numerator of the right-hand-side in 
(8) accounts for the reinforcement of the Stokes- wave am- 
plitude by the reflected acoustic wave, whereas the factor 
1 + M cos 0 ac in the denominator is attributed to the 
difference in inertial forces in the wall layer and in the 
local free-stream region. The other reason for the O(M) 
variation in the coupling coefficient as M — ► 0 is that the 
receptivity is correlated with the amplitude of the stream- 
wise velocity perturbation produced by the Stokes wave; 
the coupling coefficient in figures 1(a) and 1(b) is based 
on the mass-flux perturbation associated with the inci- 
dent wave. The amplitudes of the acoustic perturbations 
in velocity and mass flux differ by an O(M) amount that 
corresponds to the extent of density perturbation pro- 
duced by the acoustic wave. 

The coupling coefficient for both types of surface 
nonuniformities decreases as the free-stream Mach num- 
ber is increased. Throughout the range of Mach numbers 
considered, the ratio \C^ p s J. pU \/\CplJ tP u\ of the coupling 
coefficients is in the range of approximately 150 — 200 
times the ratio °f the corresponding surface- 

disturbance amplitudes. Note that receptivity that is 
induced by a given distribution of surface suction or sur- 
face temperature is generally an undesirable by-product, 
or “cost,” of using that particular type of LFC device. 
Therefore, in order to compare the costs involved in us- 
ing surface suction and surface heat transfer for LFC, it 
is necessary to adjust the ratio of amplitude pa- 

rameters in such a way that the gain in stability (i.e., the 


increase in transition Reynolds number) is the same in 
both cases. 

Because the effect on stability of a single heating or 
cooling strip is highly sensitive to the location of the strip 
(Masad and Iyer, 1993), we assume that the suction and 
the heat-transfer strip under consideration is just one out 
of many such strips that comprise the LFC system. One 
would expect that the effect on stability of this type of 
configuration will be analogous to that of continuous suc- 
tion or cooling at a suitably modified level. A compari- 
son of the results in figures 12(a) and 14(a) in Masad and 
Iyer (1993) shows that in order to achieve an identical in- 
crease in the transition Reynolds number for both types 
of LFC, the surface-temperature perturbation needs 
to be approximately 2150 times larger than the suction 
parameter The above estimate corresponds to the 

linearized limit employed in this paper, however the ratio 
«2V ^ was found to increase further as the respective 

amplitudes were increased beyond the linear range. 

For the case of = 2150, the coupling coeffi- 

cient for the heat transfer strip in figure lb will be almost 
11 to 15 times larger than the coupling coefficient for the 
receptivity produced by the suction strip. The lower ra- 
tios in the above range correspond to smaller Mach num- 
bers, whereas the higher values correspond to the tran- 
sonic range of Mach numbers. In section 3.2, we will allow 
both the frequency parameter and the orientation of the 
incident wave to vary, but keep the Mach number fixed 
at M — 0.9. The results obtained for these cases indicate 
that the ratio of the coupling coefficients for the two types 
of LFC devices varies with uj and 9 aC) but nowhere to the 
extent where it will affect the above conclusion concern- 
ing the qualitatative comparison of their respective mag- 
nitudes. Note that the estimate that = 2150 for 

equal gains in stability was based on the results obtained 
in the low-Mach-number limit, but we do not expect it to 
change significantly over the Mach-number range consid- 
ered. Of course, it needs to be seen at what values of 
and the linearized theory becomes inapprorpiate. 

3.2 Influence of Acoustic- Wave Orientation 

Now we examine the influence of the acoustic-wave 
orientation 6 ac on the receptivity produced by both types 
of surface disturbances. Toward this purpose, we consider 
the variation with 9 ac in the respective coupling coeffi- 
cients at a fixed Mach number. The location of the strip 
and the acoustic frequency parameter are the same as 
the corresponding values in figures 1(a) and 1(b). How- 
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ever, the strip geometry in (7b) is discarded in favor of 
a narrow suction and heating/cooling strip that is ap- 
proximated by the distribution = (27 t) ^ ^(^0* 

Recall from (6b) that the coupling coefficient for this ge- 
ometry is given precisely by the relevant efficiency func- 
tion ) . The reason for choosing the narrow-strip ge- 

ometry is that its Fourier spectrum is uniform across the 
wave-number range of interest and, hence, is independent 
of the wave number «ty;ib at which the geometry factor 
fti) i s evaluated. At sufficiently large Mach numbers, 
the acoustic wave number ot ac varies significantly with 
9 ac . Consequently, the choice of the geometry in (7b) 
will lead to significant variation in cv^ ; ib and hence in the 
width of the strip as the independent parameter 9 ac is 
varied across its range. This problem does not arise in 
the narrow-strip limit. Moreover, as shown by Kerschen 
and Choudhari (1985) (see, also, Choudhari 1992), a nar- 
row suction strip will result in the maximum coupling 
coefficient for to = wib, provided that the total suction 
rate across the strip is kept fixed as the strip width is 
varied. 

Because of the variation in o ac with 0 ac , the pU per- 
turbation produced by an incident acoustic wave of given 
SPL also varies significantly as 9 ac is varied across [0, 7 r]. 
In fact, for M < 1, the pU perturbation becomes zero at 
a critical angle of incidence. The critical angle is given by 
Sac — cos -1 (— M) and lies in the range (7 t/2,7t) (i.e., in 
the range of upstream-inclined acoustic waves). There- 
fore, in this section, P a c y i is kept fixed as 0 ac is varied 
(rather than fixing (pU) ac as in figures 1(a) and 1(b) of 
the previous section). 

The solid curves in figures 2(a), 2(b), and 2(c) illus- 
trate the variation in the efficiency function |A^J ;P (0 a c)| 
at M — 0.1 0.5, and 0.9, respectively. In each case, the 
efficiency function has been normalized by its value for 
9 ac — 0 (i.e., for an incident wave that propagates parallel 
to the mean-flow direction). Because of this normaliza- 
tion, the results plotted are actually independent of the 
physical quantity that is used to measure the instability- 
wave amplitude. In addition, for M — 0.1 and 0.5, the 
acoustic wave number is uniformaly small as compared 
to ai nSy [b- Therefore, in general, it will not influence the 
magnitude of the geometry factor ^^(nr^ib) to a signifi- 
cant extent, as a result of which the results in figures 2(a) 
and 2(b) also approximate the ratio of the coupling coef- 
ficients for the strip geometry defined by (7b). At each 
selected Mach number, the highest value of the efficiency- 


function magnitude corresponds to an incident wave that 
is inclined in the upstream direction. At sufficiently small 
Mach numbers (figure 2(a)), this maximum corresponds 
to an incident wave that travels in the direction that is 
nearly opposite to the direction of the mean flow. More- 
over, in all three cases, the ratio of the maximum ef- 
ficiency function magnitude to that at 9 ac = 0 is nearly 
equal to the corresponding ratio of the Stokes-wave ampli- 
tudes in the viscous sublayer, which is indicated in figures 
2(a) through 2(c) by the curve that contains symbols. In 
fact, the entire |A$.p(6^)|/|A$.p(0 ac = 0)| curve cor- 
relates with the ratio of Stokes-wave amplitudes, which 
was obtained from (8). 

The maximum of the ratio of the Stokes-wave am- 
plitudes at any given Mach number always corresponds 
to the upstream-traveling mode, and it increases rapidly 
with the free-stream Mach number. Although the maxi- 
mum of the efficiency-function curve also increases with 
Af , it increases relatively less rapidly as M approaches the 
transonic range (figure 2c). Moreover, this maximum cor- 
responds to an incidence angle of 0 ac « 145° at M = 0.9, 
and the efficiency-function magnitude decreases rapidly 
as 9 ac is increased further. Although a similar behavior 
is also noted at M ~ 0.5, it is most prominent in the 
transonic case at M = 0.9. 

The cause for the deviation of the efficiency-function 
curve from that predicted by (8) was traced to the varia- 
tions in the acoustic pressure field across the mean bound- 
ary layer, which result in a weaker acoustic-signature 
field within the viscous sublayer next to the surface. As 
remarked in section 3.1, the surface pressure coefficient 
Cp = P ac (y = 0 ) / Pq c . ; t is equal to 2 at nearly all values 
of 9 ac when the Mach number is sufficiently small. How- 
ever, this picture begins to change as the Mach number 
M becomes larger. This change is caused by a decrease 
in the acoustic wavelength (relative to the boundary-layer 
thickness) that accompanies the increase in M and leads 
to a significant refraction of the incident and the reflected 
waves across the mean boundary layer. In effect, the re- 
fraction process modifies the perfect in-phase relationship 
(i.e., a constructive interference) between the incident 
and reflected waves near the surface. This modification 
leads to a reduction in \Cp\ and also alters the nature of 
velocity fluctuations close to the surface. At a fixed value 
of the frequency parameter, the acoustic wavelength is 
relatively shorter for the upstream- travelling waves and, 
hence, the refraction by the mean boundary layer is more 
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significant for 0 ac > tt/2. As an example, consider the fre- 
quency parameter that corresponds to the lower branch 
mode at R = 955 and M = 0.9. At this frequency, 
the shortest acoustic length scale (defined as one-quarter 
of the acoustic wavelength) is approximately five times 
the displacement thickness of the mean boundary layer 
(i.e., only about 1.6 times larger than the boundary-layer 
thickness based on 99.9 percent of the local free-stream 
speed). The effect of transverse pressure variations in 
transonic boundary layers has also been noted by Ryzhov 
(1993); however, this work examined the the stability of 
boundary layers rather than their receptivity to acoustic 
disturbances. 

To check on our numerical calculations, we also com- 
puted C P as a function of 8 ac at Reynolds numbers of 
R = 1910 (i.e., R = 2.0 * 955) and R = 3820 (i.e., 
R — 4.0 * 955) for the same value of the nondimen- 
sional frequency parameter u> considered in figure 3(a) at 
R = 955. Because the mean boundary-layer profile is self- 
similar, the Cp versus 8 ac variation at all three Reynolds 
numbers was graphically indistinguishable. As a result, 
a purely inviscid calculation such as that carried out by 
Duck (1990) for M > 1 is also expected to yield the same 
result for the acoustic-field signature. 

As long as the efficiency function correlates with the 
amplitude of the Stokes shear wave, its dependence on 
0 ac is nearly identical for both types of surface nonuni- 
forrnities. (See figures 3a through 3c. for the variation of 
I p | with 0 ac .) However, if the acoustic-signature field 
becomes nontrivial because of the refraction effect at the 
higher Mach numbers, then the dependence of |A^], p | 
on Q ac tends to be different for the different type of sur- 
face nonuniformities. For example, the maximum val- 
ues of \^pij.p{8 a c\/\^ij.p{0ac = 0)| in figures 2(c.) and 
3(c) are equal to 2.2 and 3.6, respectively, for j = s and 
j = t- Similarly, in a companion study that deals with 
the receptivity induced by surface-geometry variations, 
we have found that the corresponding peak value for the 
geometry-induced receptivity is equal to approximately 
8.1. Note that the finding of a somewhat larger peak ra- 
tio in the case of temperature-induced receptivity as com- 
pared to that for the suction-induced receptivity tends to 
support our earlier remarks (section 3.1) in regard to the 
relative strength of those two receptivity mechanisms. 

The results in figures 2 and 3 can also be applied 
to predict the expected value of the efficiency function 
^pU,P 0 e > expected value of the coupling coefficient 


in the narrow-strip limit) for the case where the angle 
of acoustic incidence 0 ac is unknown, but has an equal 
probability of assuming any value in the range [0, tt] . At 
M — 0.9, the expected magnitude of A fy. p is predicted 
to be larger than |A^ ; pU (0 ac = 0)| by a factor of ap- 
proximately 1.07 for j = s and 1.54 for j = t. These 
factors are considerably smaller than the common factor 
6.15 that would have been predicted on the basis of (8). 

Thus far, we have discussed the excitation of the in- 
stability wave that corresponds to the lower branch mode 
at the location of the strip. Figures 4(a) and 4(b) de- 
pict the variation in the efficiency functions |A^. p | at 
R = 955 as the acoustic frequency parameter is varied 
while the angle 8 ac remains fixed at a given value. Note 
that a comparison of |A^. p | for j = s and j = t is 
also equivalent to the comparison of the respective cou- 
pling coefficients for an arbitrary distribution 
provided that The abscissa in fig- 
ures 4a and 4b corresponds to / = = w j and 

results have been plotted for a set of selected values for 
8 ac . For reference, the variation with / in the spatial 
growth rate — Im(ct) and the phase speed c ins = u>/Re(c*) 
is also indicated in each plot. 

For suction-induced receptivity, the magnitude of the 
efficiency function increases as the frequency parameter 
is reduced. Of course, the low-frequency modes will have 
to undergo an exponential decay before they reach their 
respective lower branch locations. Therefore, the high- 
est effective amplitude at the lower branch location will 
correspond to a frequency parameter that is somewhat 
smaller than the lower branch frequency. Alternatively, if 
one is interested in the excitation of the instability wave 
at a given frequency, then the most dangerous location 
for a narrow suction strip will be slightly upstream of 
the lower branch. For any given 0 aC) the behavior of the 
efficiency function curve |A$ ;P | is somewhat different 
from that of |A^J. p |. This difference is such that in the 
range of unstable frequencies it tends to enhance the ra- 
ti 0 |Ap(/ . p |/| Ap[/ ; p | as compared to the value of this ratio 
for the lower branch mode. On the other hand, the rapid 
rise in |A^ ;P | in the low-frequency regime will reduce the 
ratio of coupling coefficients to a small extent. Finally, 
we observe that for both types of LFC devices, the an- 
gle of acoustic incidence that leads to the most efficient 
excitation of the instability wave of a given frequency 
tends to decrease somewhat as / is increased. However, 
in general, the above orientation always corresponds to 
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an upstream-travelling acoustic wave. 

3.3 Receptivity Produced by an Array of Suc- 
tion or Heating/ Cooling Strips 

Because of the weak-surface-nonuniformity assump- 
tion utilized in this paper, the receptivity predictions for 
an individual suction or heating strip can easily be ex- 
tended to an array of strips with an arbitrary distribu- 
tion of suction velocity or surface temperature. The gen- 
eral framework to extend the Goldstein- Ruban theory for 
localized receptivity to the case of receptivity produced 
by distributed surface nonuniformities was described by 
Choudhari and Streett (1990, 1991) and in further detail 
by Choudhari (1993a,b). Basically, the efficiency func- 
tion for localized receptivity is utilized as a local influ- 
ence function in the latter case. (This idea was originally 
proposed by Tam (1981) in a related context.) The net 
receptivity is obtained by integrating over the contribu- 
tions from each individual subregion across the length of 
the region of nonuniformities. 

As shown by Zavol’skii ei al. (1983), Turnin and 
Fedorov (1983), and Choudhari ( 1993a, b), the above in- 
tegration process must include the effect of interference 
between the instability waves that are generated in differ- 
ent parts of the overall region of receptivity. The outcome 
of the interference process is determined by the geometry 
of the surface nonuniformity and by the variations in the 
wave numbers a; n5 and a ac across the region of receptiv- 
ity. In particular, when a uniform suction (or heat flux) is 
applied through an array of equidistant strips of uniform 
width, Choudhari (1992) showed that most of the recep- 
tivity is actually produced by the strips that are located 
in the vicinity of the lower branch location. The ratio 
of the effective coupling coefficient C^Xay for the array 
(which measures the effective amplitude of the generated 
instability motion near the lower branch location) to the 
local coupling coefficient c\l] al for a single strip at the 
lower branch location is given by 


n 


array 


C U) 

1 array 


K local 


F ( ' ]) {n W '\ b )\/nrD a n _, 


exp 
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(9a) 


Here a w (Rs*) = c** w 6*(Rs*) denotes the fundamental 
wave number of the periodic suction (or temperature) dis- 
tribution, and the factor D a is a measure of how rapidly 
the unsteady forcing on the right-hand side of (3) (which 


is produced by the bilinear interaction between the free- 
stream and surface disturbances) becomes detuned, or 
“desynchronized,” with respect to the phase of the insta- 
bility mode. For a flat plate boundary layer, the desyn- 
chronization factor is given by 



D<* = (a'ins-, 


lb 


, lb , 


Rs* 


(96) 


where the primes denote differentiation with respect to 
Rs * , and the subscript lb indicates evaluation at the lower 
branch location Rs* = Rs* h , as mentioned previously. In 
an accelerating or decelerating free stream, the expression 
for Da will also incorporate the additional detuning (or 
tuning) because of the variations along the streamwise di- 
rection in the physical wave number of the acoustic wave. 
Note that the ratio in (9a) is independent of the physical 
quantity in terms of which the amplitudes of the input 
and the output waves are measured, provided only that 
the same measure is adopted for both types of configu- 
rations. As a result, we have dropped the corresponding 
subscripts on in both the numerator and the denom- 
inator. 

To supplement the results of section 3.1 on the re- 
ceptivity induced by a single strip, we can examine the 
upper limit on the receptivity produced by an array of 
strips. For the strip configuration defined in (7), we 
expect that the maximum receptivity will occur when 
<y w -\b - &ins, lb - lb* In that case, the maximum 
value of the ratio in (9a) will be given by 


n 


array ; max 


lb 

fttns ; lb 


R-d 


(9c) 


where 

n D = . (9 d) 

y t7T D a 

The factor IZd is independent of 9 ac when the free stream 
is uniform. However, the ratio Rarray \ max is a func- 
tion of the acoustic-wave orientation because of the ad- 
ditional factor K a = 1 - (= 1 - for a 

two-dimensional acoustic wave) that appears on the right 
hand side of (9c). An analogous factor is also present 
in the case of receptivity to convected free-stream dis- 
turbances. However, for convected disturbances, the fac- 
tor 7 l a is independent of the disturbance structure along 
the surface-normal direction (e.g., the orientation of the 
wave-number vector if a single Fourier mode of the con- 
vected perturbation is considered). Because of an over- 
sight, the factor 7 Z a was omitted in (4a) of Choudhari 
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(1994) and, also, in the estimate for distributed receptiv- 
ity presented in Choudhari and Streett (1993), which led 
to the an overprediction of the ratio TZ array - max in those 
works. 

Not only does 7 Z a depend on the orientation of the 
incident disturbance in the case of free-stream acoustic 
waves, but can lead to either an increase or a decrease 
in the magnitude of 7£ array; max , depending on 9 ac . For 
downstream travelling waves o a c;ib is positive, which im- 
plies that \!Z a \ < 1 for those waves. On the other hand, 
o ac; ib is negative for upstream travelling waves, and can 
assume rather large values at transonic Mach numbers. 
This will lead to a substantial increase in the magnitude 
of Tlarray, max* The physical explanation for the effect 
of 7 Z a on the effective coupling coefficient for distributed 
receptivity is as follows. The length scale of distributed 
receptivity is fixed by the desynchronization factor D a . 
Therefore, it does not vary with 0 ac in the case of a 
flat-plate boundary layer. An increase or decrease in 
= rt iri 5 ,!b — o aC) ib then corresponds to a decrease or 
increase in the spacing between an adjacent pair of strips. 
The latter, in turn, implies an increase or decrease in the 
number of strips that lie within the (fixed) length scale 
of distributed receptivity and, hence, contribute to its 
cumulative effect. 

In the context of the conceptual experiment consid- 
ered in figures 1(a) and 1(b), we note that for 0 ac = 45° 
the factor 7 Z a varies from 1.0 at M = 0 to approximately 
0.87 at M = 0.99. (Because the value of 7Z a is nearly 
equal to unity in this case, the ratio 77 arrai/)ITiax could 
even be approximated by the factor 7Zo ) The magnitude 
of IZo at a lower branch Reynolds number of R = 955 
has been plotted in figure 5 as a function of M . The 
magnitude of Tip decreases from 9.3 at M — ► 0 to ap- 
proximately 7.3 as M — ► 1.0. Therefore, at 9 ac = 45°, 
the ratio |7£ array; max | will decrease from 9.3 to 6.3 across 
the range of subsonic Mach numbers. By combining this 
estimate with the results of figures 1(a) and 1(b), we can 
conclude that the effective coupling coefficients CjfJ. y 
for distributed receptivity will exhibit the same trend, as 
the Mach number is varied, as that of the respective local 
coupling coefficient for an individual strip (figures 1(a) 
and 1(b)). 

3.4 Excitation of Three-Dimensional Instabil- 
ity Waves 

The results of section 3.2 demonstrate how the de- 
crease in the acoustic wavelength as M increases can have 


a significant impact on the acoustic-signature field and on 
the magnitude of the coupling coefficients. Another con- 
sequence of the decreasing acoustic wavelength at 0(1) 
Mach numbers is that even three-dimensional instabil- 
ity waves can be excited near a two-dimensional distri- 
bution of surface nonuniformities if the free-stream dis- 
turbance has a suitable azimuthal orientation. In fig- 
ure 6(a), we have plotted the azimuthal angle <£ ins = 
arctan of the generated instability wave as 

a function of the azimuthal orientation <j> ac of the inci- 
dent acoustic wave. The Mach number chosen for this 
calculation is equal to 0.9, and the values of R and w are 
the same as the corresponding values at M = 0.9 in fig- 
ure 1(c). The polar angle of incidence 0 P C = |7t/2 — 9 ac \ is 
held fixed, and (f> ac is varied across the range 0 < <j> ac < 7 r. 
When <f) ac is an acute angle, the incident wave prop- 
agates downstream; for obtuse azimuthal orientations 
{<t>ac > tt/2), the incident wave travels upstream. Thus, 
the case = 55° and <f> ac = 180° corresponds to the 
two-dimensional wave with 9 ac — 145° in figures 1(c) and 
3(c), which leads to the highest value of the normalized 
coupling coefficient for both j — s and j = Z. 

For uj — u>ib at R — 955, instability waves with <j) ins 
less than approximately 35° can be excited near a two- 
dimensional surface nonuniformity. (See figure 6(a).) The 
maximum possible value of (f>i ns tends to increase as 9 P C 
is increased, so that the above-mentioned upper limit on 
<fii ns is realized when the incident wave propagates par- 
allel to the plate surface. At any given 1 9 P C , the same in- 
stability mode can be excited by a pair of acoustic waves 
that corresponds to two different values of <j> ac . When 
9 P C is small, the <fo n5 curve is roughly symmetric with re- 
spect to <f) ac = 7t/2. As a consequence, the wave-number 
vectors of two acoustic waves that excite the same in- 
stability mode are nearly symmetric with respect to the 
spanwise direction. As 9 P C is increased, the behavior of 
<f>ins becomes increasingly asymmetric about <f> ac = 7 t /2. 
Therefore, both incident waves, which excite an instabil- 
ity mode that corresponds to <f>i ns in the vicinity of its 
maximum, are upstream-travelling waves. The variation 
in the efficiency-function magnitudes for the two types of 
surface disturbances is shown in figures 6(b) and 6(c), 
respectively. Interestingly enough, the efficiency func- 
tion at 0 ac = 55° and 85° actually tends to zero near 
some critical azimuthal orientation that corresponds to 
an upstream-inclined acoustic wave. 
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3.5 Receptivity at Low Supersonic Speeds 

At subsonic Mach numbers, the angle 9 ac of a two- 
dimensional incident wave can be varied continuously 
from 0 to 7 r by varying its wave number a ac in the stream- 
wise direction. However, for M > 1, the acoustic disper- 
sion relation for a plane wave allows two disjointed sets of 
modes. The so-called fast mode is analogous to the acous- 
tic waves in a subsonic flow, except that the angle of in- 
cidence for a fast acoustic mode in a supersonic stream is 
constrained to the range 9 ac 6 [0, 7T-arccos ( ^)] instead 
of the complete range [0,7r] that is relevant at the sub- 
sonic Mach numbers. The excluded range of angles gets 
folded back to the range (— arccos (^) , 0) that the slow 
mode occupies. The slow mode has a critical layer within 
the mean boundary layer and is, therefore, unique to the 
supersonic range of Mach numbers. The solid curves in 
figures 7(a) and 7(b) illustrate the variation in the mag- 
nitudes of the efficiency functions |A^.p| and |A^J ; p| as 
9ac is varied for both the fast and the slow acoustic modes 
in the boundary-layer flow at M — 1.2. Again, the loca- 
tion of the LFC device is chosen to b e R = 955, and the 
frequency parameter corresponds to u> « 0.033, which 
is close to the lower branch frequency predicted by the 
quasi-parallel stability theory. Note that the effect of non- 
parallelism becomes increasingly significant (in terms of 
determining the growth rate of a two-dimensional insta- 
bility wave) as one crosses over the transonic regime into 
the supersonic range (Smith 1989, Ryzhov and Savenkov 
1990). However, the leading-order eigenfunction is al- 
ways determined by the quasiparallel theory (Fedorov and 
Khokhlov 1991); hence, the application of the present re- 
ceptivity theory is feasible in this particular case. 

Figures 7(a) and 7(b) show that the efficiency func- 
tions for the slow acoustic mode do not correlate with 
the behavior predicted by (8) (indicated by the curve with 
symbols). In fact, even the results for the fast mode show 
substantial deviation from the analytical prediction as 9 ac 
becomes large. The highest value of the normalized cou- 
pling coefficient occurs for a fast mode that is incident 
at 9 ac ^ 135° in both cases. However, its magnitude for 
temperature-induced receptivity (« 4.25) is much larger 
than that for the suction-induced receptivity (^ 1.5). The 
maximum coupling coefficient for a slow acoustic mode is 
smaller than the maximum for the fast mode in this par- 
ticular case. However, this trend is not always the case; 
because an opposite trend was encountered in an analo- 
gous set of calculations at M — 1.6. 


4. Concluding Remarks 

The Goldstein-Ruban theory was extended in this 
paper to develop theoretical predictions for the acous- 
tic receptivity produced by two-dimensional variations 
in the surface suction velocity and in the surface tem- 
perature. The receptivity mechanism that is related to 
mean-flow variations induced by these laminar flow con- 
trol (LFC) devices is likely to be stronger when surface 
heating or cooling is used for LFC (in comparison with 
the case where surface suction is used). The orientation 
of the acoustic wave can have a significant impact on 
the magnitude of the coupling coefficient in both cases. 
Moreover, at the larger subsonic Mach numbers, this in- 
fluence is rather nontrivial for upstream-oriented acoustic 
waves, which produce the highest receptivity at a given 
sound pressure level (SPL). Moreover, in the above speed 
range, the effect of acoustic-wave orientation is differ- 
ent for the different types of surface nonuniformities. Of 
course, the predictions developed in this study are lim- 
ited to small-amplitude surface nonuniformities, and fur- 
ther study is required to identify the range of validity for 
the linearized approximation, particularly in the range of 
transonic Mach numbers. 

The calculations also illustrate the influence of 
acoustic-wave orientation on the receptivity at super- 
sonic Mach numbers; however, a more detailed paramet- 
ric study is necessary to pinpoint the receptivity charac- 
teristics in this speed range, in particular because three- 
dimensional instability waves become more important at 
these Mach numbers. Such calculations will enable the as- 
sessment of the importance of receptivity near the leading 
edge and of that produced by the surface nonuniformities. 

Finally, the study demonstrates the sensitive param- 
eter dependence of acoustic receptivity in compressible 
boundary layers. This dependence will make it particu- 
larly difficult, to predict the receptivity in actual applica- 
tions and hence to predict the location of transition on a 
more rational basis. 
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Appendix 

In this Appendix, we present the inhomogeneous set 
of equations that governs the unsteady scattered field 
Qty,ac — {Cuj f acT ,aC) l ac^'-Pw ,ac} P\v,ac)' Recall that 

the superscript (j) in and Q^.ac has been dropped 
for the simplicity of notation. As stated in section 2 of 
this paper, the above set of inhomogeneous equations can 
be expressed in the following form after a Fourier trans- 
form is taken in the streamwise and spanwise directions: 

LjVs(a,<9y,/?,u;, ft(5*)Qw,<ic = N NS ;2 (Qu,,Q«c), (Al) 

where the overbars denote the Fourier transforms of the 
perturbation quantities. The operator L jvs on the left- 
hand side is the matrix form of the linear stability equa- 
tions in the coordinate system (X,Y,Z) that is introduced 
in Section 2. However, in the literature on the three- 
dimensional stability of quasi-parallel shear flows, it is 
customary to solve the relevant homogeneous equations 
in the wavefront (i.e., “Squire” coordinates) instead of 
solving them directly in the physical coordinate system. 
Hence, we will follow the same procedure here to solve the 
corresponding system of inhomogeneous equations. Thus, 
the X- and Z-momentum equations from the set (Al) 
will be replaced by the momentum equations along, and 
perpendicular to, the wavenumber vector (a, ft) in the 
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Fourier transform space. At the same time, the stream- 
wise and spanwise velocities will be replaced as depen- 
dent variables by the horizontal divergence field q+ ac = 
otCf w>ac + pW W)ac and the field q~ ac = aW w>ac - ftU w>acy 
which is related to the vertical vorticity associated with 
the unsteady scattered field. After this transformation, 
the system of equations that corresponds to (Al) may be 
rewritten in the form 


disturbance, the acoustic-signature field, and the insta- 
bility wave. Nevertheless, the magnitude of the contri- 
bution from each viscous term is already related to the 
corresponding sum of the orders of derivatives of the first- 
order perturbations Qu, and Q ac with respect to the wall- 
normal coordinate. 

The set of viscous forcing terms may then be further 
classified as 


^ J Ns(^ c ^ &Y i ft iU), Rs* )Qt£/,ae — N;v5-2(Qui)Qac)) (A 2) 

where the new array of dependent variables Q W)ac corre- 
sponds to (q+ ac , V m<ac , q~ iac ,T w<ac , P W}ae )- The operator 
L/vs represents the linear stability equations in the ro- 
tated (i.e., Squire) coordinate system and has the form 

L/vs = + B— + C , (A3) 

where the coefficient matrices A, B, and C can be ob- 
tained from the expressions given in Malik and Orszag 
(1987) and, therefore, will not be repeated in this paper. 

The forcing vector Njv 5;2 in (A2) can be expressed 
in the form 


NjYS; 2 — Nt ni ,j2 d" } (A4) 

where the vectors Nj nv; 2 and N V j sc;2 represent the bi- 
linear interaction terms that arise from the inviscid and 
viscous terms, respectively, in the Navier-Stokes equa- 
tions. The coefficients of the viscous source terms in the 
momentum and the energy equations are smaller than 
those of the corresponding inviscid terms by a factor 
of O(Rf). Hence, the error incurred in the prediction 
of the instability-wave amplitude by neglecting the vis- 
cous terms from the forcing N/^S ;2 would appear to be 
0(1/Rf), or of the same order as the error in neglect- 
ing the effects of base-flow variation (i.e., growth of the 
mean boundary layer) across the region of nonuniformity. 
However, this conjecture is incorrect because some of the 
“0(1 / Rs*)” terms involve first or second derivatives of 
the mean-flow disturbance and/or the acoustic-signature 
field with respect to Y, which can become rather large in 
certain localized regions across the boundary layer. Thus, 
the contribution from these terms (relative to that of the 
inviscid nonlinearities) is actually larger than 0(\/Rs*). 
The precise order of the contribution from each term of 
this type can only be determined after we have considered 
all asymptotic subregions that characterize the mean-flow 


Ntri 5 c;2 — + N, 


visc\2 


+ N (0 

' Vt 


sc: 2 » 


(A5) 


where the superscripts within parentheses on the right- 
hand side of (A5) indicate the sum of the orders of deriva- 
tives of and Q ac with respect to Y. A larger value 
of the superscript indicates that the overall contribution 
from that subclass of viscous terms is more dominant. 
Thus, terms that belong to N^ c;2 provide the largest 
viscous contribution, as a whole, to the amplitude of the 
generated instability wave. These terms can become com- 
parable to the inviscid forcing terms in certain localized 
regions of the boundary layer, such as the viscous layer 
close to the wall. However, as we mentioned previously, 
in the absence of a rigorous asymptotic analysis, decision 
as to whether the overall contribution from these terms 
is comparable to that of the inviscid forcing (or whether 
it simply provides the dominant correction for the latter) 
is difficult. In either case, terms from N^ c>2 provide a 
higher order correction, which is larger than 0(1/ R$*). 
On the other hand, the terms in N^ c 2 are uniformly 
of 0(1 /Rs*) with respect to the inviscid terms; there- 
fore, these terms may be safely neglected on the grounds 
stated in the previous paragraph. The only reason for 
including them in the numerical calculation would be to 
maintain consistency with the adopted framework of the 
quasi-parallel stability theory. 

The overall contribution from n£!-* , 9 and N (0 ^ 0 

ft 3 C V%SC\Z 

was small, even at Mach numbers as large as 4.5. The nu- 
merical results presented in this paper include the effects 
of both source terms; however, because the expressions 
for these terms are rather lengthy, we have omitted them 
from this appendix. Thus, only the expressions for the 
inviscid forcing terms Ni nt;; 2 and the dominant viscous 
terms are given below. The first three elements 

of each of these vectors correspond to the source terms in 
the momentum equations that correspond to q+ ac , V w>ac 
anc * iw t aci respectively; the fourth and fifth elements cor- 
respond to the forcing in the energy and continuity equa- 
tions, respectively. In deriving the expressions for these 
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source terms, we have assumed that the specific heat 
C* , the ratio of specific heats 7, and the Prandtl num- 
her <j are all constant quantities. In addition, the Stokes 
hypothesis has been assumed to be valid. The primes 
in these equations denote derivatives with respect to the 
wall-normal coordinate Y. Because of the lengthy na- 
ture of the expressions involved, the latter have also been 
checked with the symbolic manipulation package Mathe- 
matica. Expressions for N^2c;2 and are available 

from the author upon request. 

Nj„„; 2 (l) = ^-{-iup w {aUae + pW ac ) 

Vo L 

-4 -—[i(°tacUw 4 fiac.W w )(otUac 4 PW ac) 

f 0 

+V w (*U' ae + f3W^ c ) 

+i(c* w Uac 4 fix V W a c)( a Pxv 4 fiW w ) 

+V ac (aUi +0K)} 

+Pui [*(«ac^0 + /?aeWo)(<*l^ae + 0W a c) 

+v ac ( aUl + pWD] 

+Pac [*(<*«; t^O + Pw Wo){aU w + pw w ) 

+V w {aUl+0Wo)]} 

3 Rs* ( 

Ninv;2(2) = - — {-iupwV ac 

4 vo L 

-f— [i(a ac U w 4 fiacW w )V ac T Vtv Vac 
Jo 

4 i(ot w Uac T PwWac)Vw 4 KicK^] 

+ipw(<XacUo + fiacWo)Vac 
+ip ac (a w U Q + 0 w W o )V w \ 

Nin«;2(3) = — \-iuip w (-PU ac + OtW ac ) 

+ ^r[i(a ac U w + PacW w )(-pU ac + aW ac ) 

Jo 

+V w (-pU' ac +aW' ac ) 

+i(a w U a c + PwW ac )(-pU w + aVVu,) 


+V ae (-pUl+aWi)] 

+Pw [t(«ac^0 + PacWo)( — pU a c + OtW a c) 

+V ae {-/3Ul + aWD] 

+Pac [*(«tu Uq + Pw Wq)(—0U w + 

+ V w (-pU^ + aWti]} 

N>ni<;2(4) = / — iu>p w T ac 

Vo K 

4;pr \K a ^cPw 4 fiacW w )T ac 4 Y w T ac 
Jo 

+ i(a w Uac 4 A vW ac )T w 4 VqcT^] 

~\~Pw [*(^acf^0 4 fiacWo)T ac 4 V^ c Tq] 

4 pac \i( a ™Uo 4 ftw Wo)T w 4 Ku^o] 
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Figure 1(b). Influence of Mach number on magnitude of \Cp U;p u\/ € w 
T stgn = 311 K. R = y/RT t = 955. w = a> lb . /? ins = 0. and 





Figure 2(a). Influence of acoust ic-vave orientation on |A^. p | with 


M = 0.1. T lan = 311 K, R = x/ReZ = 955, 


cj ^ CJ, 


and 






Figure 2(b) . 


Influence of acoustic-wave orientation on |A^]. P | with 
M = 0.5, r rfyTl = 311 K, R = = 955 - ^ * w lb , and 
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Figure 3(b). 


Influence of acoustic-wave orientation on |A^y. p | with 
M = 0.5, = 311 K . R = \JRe 2 = 955, o w lb , and 
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Figure 3(c) . 


Influence of acoustic-wave orientation on |A^y.p| with 
M = 0.9. T stgn = 311 K, R = \/Kel = 955, cj w lb , and 

Pin = o.o. 
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f=cj'v'/U*J 


Fig. 4(a). Variation of |A£y. p | with / at fixed strip location 
with = 311 K. R = JRel = 955, and /?* = 0. 





/=cjV/LC 2 


Fig. 4(b). Variation of |A^. p | with / at fixed strip location 
with T tm = 311 K, R = yJRe'i = 955. and 0^ = 0. 






Figure 6(a) . Excitation of three-dimensional instability waves 

near two-dimensional nonuniformity : 0^ (deg) is 

plotted against 0^ (deg) with U = 0.9, = 311 K, 

R = \fKe~ 2 - 955, and / ; u/R S ‘ = 18E-06 . 
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Figure 6(b). Excitation of three-dimensional instability waves 
near tvo-dimensional suction strip: |A^y. p | is 

plotted against 0^ (deg) with id = 0.9, T sm = 311 
R = \JRe 2 = 955, and f = cj/R s - = 18E-06 . 
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Figure 6(c). Excitation of three-dimensional instability waves 

near two-dimensional heat-transfer strip: |Apy.j>| is 

plotted against 0^ (deg) with U = 0.9. = 311 K, 

R = y/Re~ t = 955, and / = = 18E-06 . 
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Figure 7(a) . 


Influence of acoustic-wave orientation on 
U = 1.2, T stgn = 311 K. R = = 955, u * cj lb . 


180 


| with 
and 
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Figure 7(b) . 


Influence of acoustic-wave orientation on with 

M = 1.2, T stgn = 311 K, R = = 955, u =« w lb , and 
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